UNITED STATES SPECIFICATION 



TO ALL WHOM IT MAY CONCERN: 

BE IT KNOWN THAT WE, JORG KLIEWER, citizen of Germany, 
residing at Nettelbeckstrasse 4, 24105 Kiel, Germany and 
Georg WEISS, citizen of Germany, Haupstrasse 15, 85110 
Kipfenberg, Germany, have invented certain new and useful 
improvements in a 

METHOD FOR THE LINEARIZATION OF FMCW RADAR DEVICES 
of which the following is a specification. 
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CROSS REFERENCE TO RELATED APPLICATIONS 
This application is based upon German Application 
Serial No. 10315012.9 Filed on April 2, 2003 wherein 
priority is claimed under 35 U.S.C. 119. 

BACKGROUND OF THE INVENTION 

The invention relates to a method for the linearization 
of a frequency modulated continuous wave (FMCW) radar 
device . 

With FMCW radar devices, phase errors occur in the 
reception signal, as the result of a not strictly linear 
increase in a transmission frequency of an FMCW radar 
device. This results in a loss of accuracy in the 
evaluation of the signal. 

In this case, various methods are known for linearizing 
an FMCW radar device. For example, one known method 
discloses a method wherein a delayed transmission signal is 
generated via an optical delay line having a defined running 
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time. In this case, the delayed transmission signal is 
mixed with the non-delayed transmission signal. Thus, it is 
possible to determine the deviation from the ideal linear 
characteristics. In particular, the disadvantages of this 
method are the great amount of construction effort, expense 
and the age-related clouding of the optical components. 

In addition, it is also known to generate a signal 
function that can be predetermined, such as a linear 
characteristic via a digital direct synthesis (DDS) system. 
However, a disadvantage of this system is that there is a 
large construction effort and expense associated with this 
cost intensive structure. 

An FMCW radar system having linear frequency modulation 
is known from German Patent DE 4104792 Al . In this system, 
there is an oscillator that generates a ramp-shaped 
modulated signal having a modulated oscillation frequency, 
which is emitted by an antenna. A monitoring and control 
unit determines the deviation of the measured modulated 
oscillator frequency of a signal reflected from an object 
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via a Hilbert transformation using a linearized reference 
value. In this case, a phase deviation can be determined 
from the deviation that has been determined. 

From this deviation, a correction signal is generated, 
which is passed to the oscillator, whereby the oscillation 
frequency is modified in accordance with the predetermined 
reference values. A disadvantage of this system is that 
this involves a great construction effort, which is cost 
intensive, since additional components are required. 

SUMMARY OF THE INVENTION 

The invention relates to a method with which an FMCW 
radar device can be linearized, without a great construction 
effort or expense. 

In this case, a correction phase term is calculated on 
the receiver side for compensation of the phase error in the 
reception signal q(t). In contrast with known methods, no 
costly components are required in the transmission branch of 
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the FMCW radar device when using this method and this thus 
results in reducing costs. 

This method, which relates to the invention should 
include the following steps: 

selecting a number L of consecutive ramp-shaped 
reception sequences qk(n) of the reception signal, which 
number can be predetermined, with K=1,...,L, the phases 
arg{q k (n) } of which can be represented as a polynomial of the 
N th order for the time index n, with the polynomial 
coefficient m*, with {=1,..., N, 

transforming a spectrum range Q(e jfi )of the selected 
reception sequences q(n) into a basic band reception 
sequence q k {n) with k=0,...L-l> 

iteratively calculating a correction phase term for 
partial compensation of non-linear frequency components in 
the basic band reception sequences q k ( n ) by means of 



4 



calculating polynomial coefficients in\[l of the individual 
basic band reception sequences via estimating methods 

wherein are the sequences that have already been 

iteratively phase-corrected, and whereby the iteration is 

stopped once the parameter change between two consecutive 
iterations, which can be predetermined, remains below a 
threshold c, which can be predetermined. 

In another embodiment of the invention, the scanning 
cycle T A of the ram signals is reduced, in a method step 
after the basic band transformation. For this purpose, a 
prior Antilias low pass filtering of the ramp signals q k (n) 
is performed. In addition, it is beneficial if the factor 
K, by which the scanning cycle T A is reduced, lies between 
K=30 and K= 60. 

In contrast, with an alternative linearization of the 
frequency modulation of the radar transmission signal, the 
linearization is performed on the receiver side, via a 
digital processing of the reception signal. Thus, no 
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additional components are required for this method. The 
calculations required for this method can be performed by 
the components for signal evaluation that are already 
present. This results in additional advantages with regard 
to the production costs. 

Another advantage of the invention is that there can be 
in an inexpensive estimation method that can be easily 
implemented, since the reception signal is processed in a 
basic band. This signal results in a reduction of the 
scanning cycle T A by a factor K between 30 and 60. In 
addition another advantage of this method is that there can 
be a calculation of the polynomial coefficients via a robust 
estimation method. This is even in the case of channels 
subject to severe interference. Thus, there can be a 
relatively good compensation of phase errors with signal to 
noise ratios (SNR) of up to -20 dB using this method. Thus, 
the method can be used for not only the civil sector, but 
also in the military sector or for distance measurements. 



BRIEF DESCRIPTION OF THE DRAWINGS 



Other objects and features of the present invention 
will become apparent from the following detailed description 
considered in connection with the accompanying drawings, 
which discloses at least one embodiment of the present 
invention. It should be understood, however, that the 
drawings are designed for the purpose of illustration only 
and not as a definition of the limits of the invention. 

In the drawings, wherein similar reference characters 
denote similar elements throughout the several views: 

FIG. 1 is an exemplary idealized ramp-shaped FMCW 
transmission signal ; 

FIG. 2 is an exemplary representation of the refection 
of the radar signal on a point shaped object at distance R; 

FIG. 3 is a schematic diagram of the method that 
relates to the invention; 
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FIG. 4 is a representation of a basic band shift of a 
selected spectrum range in the reception signal; 

FIG. 5a is a representation of a spectrum of a cycle- 
reduced basic band signal before and after the band pass 
filtering; and 

FIG. 5b is a representation of an amount spectrum of a 
radar reception signal before and after the phase 
correction . 



Turning now in detail to the drawings , FIG. 1 shows the 
starting point of this method is the transmission signal of 
an FMCW radar that is in the general polynomial approach of 
the N th order for a transmission ramp as: 



with 0=2,..., N represent the polynomial parameters, yo 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 




(1) 



With fe[0,r r ]. With this design, the coefficients m c 



represents the initial phase, and T r represents the time 
duration of the transmission ramp. The parameter mi or slope 
that is assumed to be given, indicates the linear component 
mi=B/t, which is exclusively present in an ideal case, 
wherein B describes the frequency stroke of the FMCW radar 
system. Thus in the following formula, a single reflecting 
point shaped object is assumed. Thus, the corresponding 
parameters are represented in FIG. 1, using the example of 
an idealized transmission signal. Thus with Equation (1) 
the reception signal in an ideal, noise free case is then: 

y(t) = exp{j27i^fj + jt 2 +—t 3 +— t 4 +...+y Q -j^j} (2) 

Since the determination of relative distances between 
objects that are moving relative to the radar sensor at 
approximately the same speed represents a possible use of 
the method described, it can be assumed that s(t) = 2R in 
the following (Fig. 2). 
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For the momentary wavelength Mt) , the following 
results at the speed of light c: 




/o+2>/] with /(/) = 



1 flfarg{x(Q} 
In dt 



It is advantageous if the reception signal y(t) is mixed 
into a lower frequency position. The mixed reception signal 
q e (t) results as follows, from Equations (1) and (2): 



If q e (t) describes a harmonic exponential function, i.e. 
if m* * 0 with { = 2, . . . , N, then the ideal case is present, 
and the transmission signal x(t) possesses a purely linear 
increase in the transmission frequency. 

In order to indicate the maximal deviation from this 
ideal transmission frequency progression, a maximal relative 
linear error is defined as 



q e (t) = y(t) • x (0 = exp - j — fo + ±,m e f 

c V f>=\ J 




(3) 
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with fi d (t ) =fo+mit for the ideal case. The symbol 0 
therefore describes a maximum deviation of the momentary 
frequency progression for f 0 = 0 from the linear function 

fidiOlf-Q- Taking into consideration the noise term r(t) that 

is always present in real applications, the following is 
obtained from Equation 3, with q(t) = q e (t)+r(t) 

q(ty= expj- J^{fo + Z + KO (4) 

If the complex-value noise signal r(t) is assumed to be 
white and gauss-distributed in Equation 4, for the sake of 
simplicity, the noise influence for a significantly large 
signal-to-noise distance can also be considered as phase 

noise e jr{t) , whereby ^(0 is real-value, white, and gauss- 
distributed. The phase signal, which is made time-discrete 
with the scanning cycle T A , results from 
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arg{q(nT A )} = (f>(nT A ) 



4ttR 



N \ 



C 



fo + Hrn e {nT A Y + r(nT A ) 



(5) 



Proceeding from the known (received) phase scanning 
values <t>(nT A ), the method according to the invention 
implements a robust estimate of the polynomial parameters m c , 
with 0 = 2, . .., N. In an advantageous embodiment of the 
invention, a robust estimate of the unknown object distance 
R is additionally possible. 

Fig. 3 shows a block schematic in the form of a 
sequence schematic of the method according to the invention. 
The initialization 1 takes place, according to the 
invention, by means of a basic band shift BV of L 
consecutive ramps of the sequence q(n) = q(nT A ) that results 

from Equation (4), by means of scanning. The shift of a 



into the basic band is shown as an example in Fig. 4. The 
cut-off frequencies Qi and Q u are defined, for example, by 
way of the minimum and maximum distances R min and R ma x/ which 



spectrum range, which can be predetermined, 
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are determined by the purpose of use of the FMCW radar 
system: 

4tv 4k 1 

c c 2 

The ramp signals qjc(n) that are formed in this manner 
possess a band-pass spectrum. It is advantageous that after 
the basic band shift and the subsequent Antialias low-pass 
filtering TF, the scanning cycle can be reduced by a factor 
K. In this connection, the factor K particularly lies 
between K = 30 and K = 60. This results in advantages with 
regard to a reduction in effort/expense of all further • 
partial operations . 

The next step or progression of the block diagram in 
Fig. 3 shows the sequence of the iterations for the 
calculation of a correction phase term for partial 
compensation of non-linear frequency components in the basic 
band reception sequences q^im) with m as a time index after 

cycle reduction. In the further process steps, a signal 
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sequence whose scanning cycle T A was reduced is assumed with 

q { p{m) . Of course, these process steps can also be carried 
out with a sequence ^(ri), whose scanning cycle T A was not 
reduced. 

Within the loop that runs through in each iteration, it 
is advantageous if band-pass filtering of the sequences 

q^ijri) is located at the beginning, which is advantageously 

carried out as window formation in the spectrum range, 
whereby all the spectral coefficients outside the window are 
set to zero. It is advantageous if a Hamming window is used 
as the window function. 

In the case of several reflecting objects, the 
reception signal spectrum of a single object is selected by 
means of the band-pass filtering. All of the other, 
undesired spectra are suppressed. In addition, suppression 
of those noise-like spectrum components that do not belong 
to the working component takes place. In this way, the SNR 
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can be increased after the band-pass filtering, even in the 
case of a single reflecting object. 

To calculate the band-pass filtering, it is 
advantageous to perform a fast Fourier transformation FFT of 
the sequences q { p(m) . For this purpose, it is advantageous 

if the individual discrete Fourier transformed Q^iju) of the 

basic band reception sequences q { p(m) are calculated, whereby 

Q¥ ) Qi) = FFTffi\m)} for k = 1,..., L. The point count M of the 

FFT typically lies at up to 1000. In the case of a 
reduction of the scanning cycle T A that has advantageously 
been performed, and prior Antialias low-pass filtering, the 
FFT can be carried out with a longer point count M of 64 to 
256. 

In a next iteration step, filtered basic band reception 
sequences q^\m) are calculated by means of a band-pass 

filter according to &<'>(//) = ^)&' y (j*) > with & (,) (//) = FFT{q^{m)) 
, whereby w(p) is a spectrum window that can be 
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predetermined, and //£[//„,//,] indicates the range of the 

spectrum window having a window center point |a max that can be 

predetermined, whereby, with a lower limit \x u that can be 

predetermined, and an upper limit |i i that can be 
predetermined. In order to determine the position of the 
center point |i m ax, the FFT amount spectra are averaged 

over the number L, and the maximum is sought from the 
averaged spectrum, according to the following formula: 

For example, for a single reflecting object, jlx i = 0 and 
ji u = M-l. For several non-overlapping objects, the limits |i i 
and |x u must be selected so that in each instance, only the 
spectral component that corresponds to the object of 
interest is detected. 
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Subsequently, it is advantageous if the calculation of 
the individual inverse discrete Fourier transformed q^{m) of 

the filtered FFT spectra Q^i/u) takes place, whereby 
ql i) (m) = IFFT{Q k {i) (M)} for k = lf L " 

In a subsequent iteration step, the distances Rk {l) are 
estimated by means of a maximum likelihood estimating 
method. The polynomial coefficients Rk 0) can be estimated 

from the estimated distances m a . Subsequently, the 

distances R k are averaged over L reception sequences q k (n) . 

In a final iteration step, a reception sequence q ( k i+]) (m) with 

the averaged estimated polynomial coefficients m i is 

calculated as the starting point for a next iteration loop, 
according to 
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The maximum likelihood parameter estimate takes place by 
using the linear estimator 



v = [A T A]' l A T b (6) 

with the vectors v = [Rf,Jm x Jm 2 ,..^> 6 = [<D(0) 5 O(r 4 ) 5 O(2rj ? ...] 7 and 

a matrix A having the dimension L q x (N+l) , whereby L q 
represents the length of a sequence q { k °(m) . The 

determination of the phase values <|>(nT A ) in the vector b 
takes place by means of a phase unwrapping method. 

A lesser linearity error in the cycle-reduced basic 
band signal sequence q { k n (m) is successively achieved by means 

of the phase correction, and this results in a lesser 
estimate variance, in each instance. If the condition 

|^ ( '" !) - ^ (/) | < £r has been fulfilled after i=l iterations, final 
estimated values R and m t can be calculated from the 

individual intermediate results R and '■ , which have been 
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averaged over all L reception sequences in the i iteration, 
in each instance, according to: 

R = R U) and m, = ^E* (0 «{° ,^=2,..N 

In the following, it will be shown that it is possible, 
using the method according to the invention, to reconstruct 
the spectral resolution, which has been reduced as a result 
of the linearity error in the FMCW radar transmission 
signal, of the scanned, mixed reception signal 'q(n). This 
will be shown in exemplary manner below, using the example 
of two objects at different distances. 

The starting point is a quadratic transmission signal 
sequence (chirp) with a linearity error of 0=5%, whereby 
the simulation parameters are selected to be N = 2, m 2 = 
5xl0 12 , SNR = -18 dB, K = 40, and M = 64 . For the other 
system parameters, the following table applies: 
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\/2 1 lip 


Minimum possible distance 


R min = 900 m ! 


Maximum possible distance 


R max = 1,100 m 


Ramp increase time 


T r = 1 ms 


Basic frequency 


f 0 = 35 GHz 


Frequency stroke 


B = 100 MHz 


Scanning frequency 


f A = 2 MHz 


Spectrum averaging L 


L = 20 


Spectrum averaging, total 


25 x L = 500 


Spectrum window w(ji) 


Hamming window 



The objects are located at the distances R 0 = 900 m and Ri = 
1,100 m. Fig. 5 shows the amount spectra 

the cycle-reduced signal q^\m) for a fixed k = 0, L-l 

during the first iteration before (top representation in 
Fig. 5a) and after (bottom representation in Fig. 5b) band- 
pass filtering. The partial spectra, which can be easily 
separated, are clearly evident, whereby the bottom partial 
spectrum, which belongs to the object at the distance Ri, is 
selected in the band-pass filtering. 



The top figure in Fig. 5b shows the range of the amount 



spectrum 



of the radar reception signal that is of 



interest, which is established by the minimum and maximum 
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distance (see table) . As the result of the non-linearity in 
the transmission signal, the frequency lines that correspond 
to the two objects are smeared and can hardly be recognized. 
Using the parameters estimated for the object at the 
distance Ri, using the method according to the invention, the 
inherently non-linear component of this signal is now 
compensated, as- shown in the bottom part of Fig. 5b. 



A comparison of the figures shows the clear increase in the 
object resolution. The small triangles in the bottom part 
of Fig. 5b characterize the positions of the frequency lines 

in the amount spectrum j^^ 7 ") , which would generate a 
transmission signal having a linear frequency increase. 



Accordingly, while at least one embodiment of the 
present invention have been shown and described, it is to be 
understood that many changes and modifications may be made 
thereunto without departing from the spirit and scope of the 
invention as defined in the appended claims. 
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